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1 Introduction 

We present a study of the intermediate regime between ultra-relativistic and 
nonrelativistic flow for gamma-ray burst afterglows. The hydrodynamics of 
spherically symmetric blast waves is numerically calculated using the AMR- 
VAC adaptive mesh refinement code. Spectra and light curves are calculated 
using a separate radiation code that, for the first time, links a parametrisa- 
tion of the microphysics of shock acceleration, synchrotron self-absorption 
and electron cooling to a high-performance h ydrodynamics simulation . 
All results are presented in more detail in IVan Eerten et al. I (|2nn9l ). 



2 Dynamics 

We used typical physical parameters: explosion energy 10^^ erg, circumburst 
number density 1 cm~^, assuming a homogeneous environment, shocked 
magnetic field energy density a fraction 0.01 of thermal energy density, 
shocked electron energy density a faction 0.1 of thermal energy density. We 
trace the evolution of the magnetic field and gradually decrease the frac- 
tion of the number density for accelerated particles at the shock from 1 
downto 0.1. 

The blast wave slows down slower than anticipated from equating the 
explosion energy to the total swept-up rest mass energy, which in this case 
would lead to a transition time tNR ^ 450 days. We find ^ 1290 days in 
the emission time frame. The analytical est imate for blast wave velocity x 
Lorentz factor (Pj) from Huang et al. ( 19991 ) overestimates the late time Pj 



by a factor 4/3. 

The simulation uses an advanced equation of state that resulting in an 
effective adiabatic index Fad gradually evolving from 4/3 to 5/3 from rel- 
ativistic to nonrelativistic. We find that the lab frame density directly 
behind the shock D, divided by po7^ (where po the circumburst den- 
sity) remains numerically close to four throughout the entire simulation: 

D/pol^ = (Tad + l/7)/(rad - 1) - 4. 
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Figure 1: Left: Spectra at different observer times. The smooth curves show simulated 
spectra at different observer times: 1, 10, 100, 1,000 and 10,000 days, with later observed 
spectra having lower flux in the high frequency range. For comparison we have included 
predicted slopes at the different power law regimes after 1 day, for both ^jv = 1 (solid line) 
and = 0.1 (dashed hne). Right: Simulated light curve at 200 MHz for GRB030329, 
top curve for spherical explosion and bottom curve for hard-edged jet with opening angle 
22 degrees. We have drawn the following slopes from left to right: 1/2, 5/4 and 11/10. 
LOFAR sensitivity for 25 core and 25 remote stations after four hours of integration time 
is 0.273 mjy and indicated by the horizontal line. The contribution from the counterjet 
is also shown separately. 

3 Radiation and GRB030329 

Our code successfully reproduces the synchrotron spectrum, including self- 
absorption and electron cooling (see fig. [1]). The transition to the non- 
relativistic regime occurs around 1000 days in observer time. Only the 
relativistic slopes of the light curve match the analytically expected values. 
The changing Fad leads to a late-time steepening, but the light curve is then 
more strongly influenced by the change in ^n, leading to a less steep decay. 

Simulating GRB030329 using explosion parameters derived from analyt- 
ical model fits reveals a strong discrepancy between simplified model and 
more accurate simulation, with fiux differences up to an order of magnitude. 
Simulated late time radio light curves using hard-edged jets indicate that 
the counterjet should eventually become visible, and show a doubly peaked 
curve due to the combination of strong self- absorption and the jet break. 
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